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Part  1.  A  summary  of  the  Ph.  D.  thesis  of  Jinghua  Shi  entitled 


QUANnnCATION  OF  EARTHQUAKES  BASED  ON  IMPROVED  ESTIMATES  OF 

ENERGY  USING  REGIONAL  WAVES 


INVESTIGATIONS  UNDERTAKEN 

We  report  our  studies  of  source  and  wave  propagation  parameters  for  regional  earth¬ 
quakes,  especially  the  precision  in  measurement  of  earthquake  strength,  using  seismic  energy  as 
the  underlying  property.  The  study  is  accomplished  by  analyzing  the  strongest  regional  seismic 
waves  at  distances  in  the  range  20  -  1000  km,  using  frequencies  across  a  wide  band.  Because  of 
the  demonstrable  stability  of  root-mean-square  (RMS)  Lg  amplitude  and  the  relationship  between 
RMS  Lg  and  seismic  energy,  the  study  improves  the  accuracy  of  the  traditional  steps  of  correcting 
measured  regional  wave  amplitudes  for  propagation  effects,  such  as  geometrical  spreading  and 
attenuation  (1/0  as  a  function  of  frequency,  and  on  the  scaling  of  regional  earthquakes  in  north¬ 
eastern  United  States,  such  as  static  stress  drop  and  seismic  energy. 

RESULTS 

We  have  carried  out  three  projects: 

•  The  method  of  spectral  fitting  a  oo-square  source  model  to  the  Lg  displacement  spectrum  is 
applied  to  make  a  more  reliable  source  comer  frequency  and  stress  drop  estimates. 

•  RMSLg  is  defined  from  vertical  component  Lg  waves,  the  quantity  which  has  its  underlying 
physical  significance  in  terms  of  seisnuc  energy  and  shows  excellent  consistency  and  stability 
for  single  station  measurement. 

•  Based  on  the  close  relation  between  RMS  Lg  and  the  radiated  seismic  energy,  the  empirical 
scaling  of  the  total  seismic  energy  (calculated  from  the  source  parameters  with  an  assumed 
rapture  model)  with  the  network  averaged  RMS  Lg  estimates  is  obtained  both  observationally 
and  theoretically. 

The  Comer  Frequencies  and  Stress  Drops  of  Intraplate  Earthquakes  in  the  Northeastern  United 
States 

This  study  presents  the  estimation  of  stress  drops  for  small  to  ntiddle  size  intraplate  earth¬ 
quakes  in  the  northeastern  United  States.  The  vertical  component  Sg  and  Lg  waves  of  49  earth¬ 
quakes  were  analyzed.  From  our  study,  the  source  comer  frequencies  estimated  by  fitting  the  Lg 
displacement  spectrum  with  the  assumed  co-square  source  model  are  more  consistent  with  the  cor¬ 
ner  frequencies  measured  from  empirical  Green’s  function  deconvolution  method  than  those  esti- 


mated  from  the  cross  of  horizontal  low-frequency  spectral  asymptote  and  a  line  indicating  the  (0‘^ 
decay  after  comer  frequency  (about  35%  systematically  smaller  than  the  comer  frequencies  mea¬ 
sured  from  empirical  Green’s  function  method).  The  source  comer  frequencies  we  estimated 
turned  to  be  most  appropriate  for  the  small  to  middle  size  earthquakes. 

The  estimates  of  the  stress  drop  for  most  of  the  events  fall  between  0.1  and  200  bars  (Fig¬ 
ure  1).  The  stress  drops  for  events  with  the  seismic  moment  less  than  about  2  x  10^”  dyne-cm 
increase  with  increasing  moment,  then  become  less  dependent  on  moment  for  the  larger  events. 
The  critical  value  of  2  x  10^°  dyne-cm  seismic  moment  corresponds  to  nearly  about  magnitude  3. 
A  dependence  of  stress  drop  on  moment  for  smaller  events  suggests  a  breakdown  in  self-similar- 
ity  below  a  threshold  moment. 

The  Stability  of  RMS  Lg  Values  in  Estimating  the  Size  of  Regional  Earthquakes  in  New  York 
State  and  Adjacent  Areas 

Here,  we  investigate  the  property  of  Lg  waves  in  terms  of  RMS  Lg  (root  mean  square 
value  of  Lg  wave  ground  velocity  amplitude  within  a  certain  group  velocity  range  and  scaled  to 
100  kilometer  reference  epicentral  distance),  especially  the  stability  of  RMS  Lg  values  for 
natural  earthquakes,  in  order  to  develop  a  magnitude  scale  that  improves  upon  current 
procedures  based  on  coda,  or  body  waves,  or  maximum  amplitudes. 

From  52  regional  earthquakes  in  the  Northeastern  United  States,  the  RMS  Lg  values  we 
obtained  for  epicentral  distances  (A)  up  to  900  kilometers  decrease  with  A  as  about  after 
correcting  for  anelastic  frequency  dependent  attenuation.  The  measured  RMS  Lg  values  from  dif¬ 
ferent  stations  are  quite  consistent  to  each  other,  indicating  excellent  stability,  as  shown  in  Figure 
2.  RMS  Lg  values  from  single  station  measurement  are  consistent  between  different  stations  and 
are  quite  stable  for  regional  earthquakes.  The  network  averaged  Log\Q(RMS  Lg)  against  seismic 
moment,  LogiffiM^,  result  in, 

LogyQ^RMS  Lg)  [nm/s)  =  Q.^lLogy^iM^  [dyne-cm)  -12.70 
with  orthogonal  standard  deviation  of  0.12  Log^fMo)  unit. 

Determination  of  Seismic  Energy  from  Lg  Waves 

Lg  waves,  which  can  be  treated  as  wide  angle  reflection  of  shear  waves  trapped  within  the 
crust,  are  the  dominant  feature  of  seismograms  for  regional  earthquakes,  and  become  very  useful 
for  regional  earthquake  analysis,  specifically  for  the  estimation  of  radiated  seismic  energy.  In  this 
study,  we  make  the  estimates  of  seismic  energy  from  regional  Lg  waves,  using  broadband  fre¬ 
quency  data  of  regional  earthquakes  in  a  region  where  the  attenuation  information  is  well  deter¬ 
mined.  We  have  demonstrated  the  seismic  energy  radiation  for  different  classes  of  earthquakes 
using  Lg  waves  in  terms  of  RMS  Lg,  and  obtained  the  scaling  of  the  total  seismic  energy  of  earth¬ 
quakes  (using  an  assumed  source  model)  with  seismic  energy  measured  from  Lg  waves  (RMS  Lg) 
both  synthetically  and  observationally  for  regional  earthquakes  in  the  northeastern  United  States. 

The  relation  between  the  total  seismic  energy  calculated  from  the  earthquake  source 
(assuming  Boatwright’s  (1980)  Deceleration  model)  and  the  network  averaged  RMS  Lg  value  can 
be  expressed  as  (Figure  3), 

For  52  regional  earthquakes: 


[ergs]  =  2,0^log^Q(RMS  Lg)  [nm/s]  +8.91 


a  =  0.15 


From  the  synthetics: 
for  vertical  strike-slip, 

iogiQiE^)  [ergs]  =  llog^^^RMS  Lg)  [nm/s]  +9.014  STD  =  0.03  ; 
for  30°  thrust  fault, 

logiQiE^)  [ergs]  =  llog^^^RMS  Lg)  [nm/s]  +8.718  STD  =  0.03  ; 
for  60°  normal  fault, 

logy^iE^)  [erg^]  =  llog^^iRMS  Lg)  [nm/s]  +8.715  STD  =  0.03  . 

The  seismic  data  upon  which  our  results  are  based  are  available  in  AH  format.  The  contact  per¬ 
son  is  Jinghua  Shi  (phone:  914-365-8460;  Email:  jinghua@ldeo.columbia.edu). 

NON  TECHNICAL  SUMMARY 

In  this  project,  we  studied  the  use  of  Lg  waves  in  terms  of  the  precise  scaling  of  regional 
earthquake  strength,  such  as  the  stress  drop  and  seismic  energy  release.  We  defined  a  quantity, 
RMS  Lg  (root-mean-square  of  Lg  amplitude  scaled  to  100  kilometer  reference  distance)  wide  fre¬ 
quency  integration  of  vertical  component  Lg  signals  within  certrdn  group  velocity  window.  We 
proved  that,  RMS  Lg  has  a  physical  significance  in  terms  of  seismic  energy,  shows  excellent  con¬ 
sistency  and  stability  even  for  single  station  measurement.  Based  on  the  close  relation  between 
RMS  Lg  and  the  radiated  seismic  energy,  the  absolute  scaling  is  obtained  between  the  total  seis¬ 
mic  energy  from  earthquake  source  and  the  network  averaged  RMS  Lg  estimates. 
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Figure  1.  The  relation  between  the  static  stress  drop  and  the  seismic  moment  with  the  error 
bar  (the  mean  of  error  in  the  event  network  average).  The  filled  circles  are  for  the  events  for 
which  only  short-period  data  are  available.  The  empty  circles  are  for  the  events  that  have 
broadband  data.  The  triangles  are  for  the  events  whose  source  comer  frequencies  were 
estimated  from  empirical  Green’s  function  method. 
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Figure  3.  The  regression  of  seismic  energies  calculated  from  the  source  parameters  with 
the  network  averaged  RMS  Lg  values  for  52  regional  earthquakes  in  northeastern  United 
States  is  plotted.  We  assumed  the  source  as  Boatwright’s  Deceleration  model  in  seismic 
energy  calculation.  The  network  averaged  RMS  Lg  values  are  measiured  from  vertical 
component  Sg  and  Lg  signals  within  a  group  velocity  window  from  3.6  -  3.0  km/s. 


Part  2.  An  application  of  methods  developed  by  the  thesis,  described  in  Part  1,  to  estimate  the 
spatial  variability  of  seismic  wave  attenuation  in  Western  China. 


Attenuation  Variation  in  Xinjiang,  China 

Jinghua  Shi,  Paul  G.  Richards,  and  Liping  Gao 
Lamont-Doherty  Earth  Observatory  of  Columbia  University 
Palisades,  NY  10964 

ABSTRACT 

We  have  applied  the  empirical  Green’s  function  method  to  event  pairs  recorded  by 
the  single  station  WMQ  (Urumqi,  China)  from  1988  -  1989  to  determine  the  source  comer 
frequencies  of  the  28  larger  events  in  these  pairs,  and  have  estimated  the  crustal  averaged 
attenuation  quality  factor  Q  around  station  WMQ.  28  path-averaged  quality  factor  Q  val¬ 
ues  were  obtained  from  epicentral  distances  ranging  between  90  to  260  kilometers.  We 
found  that:  the  Q  values  for  the  paths  crossing  the  Tien  Shan  mountains  are  the  lowest, 
about  500  -  600;  the  Q  values  for  the  paths  along  the  Tien  Shan  mountains  are  relatively 
higher,  about  700;  and  the  Q  values  for  the  paths  through  the  basin  to  the  southeast  of 
WMQ  are  the  highest  in  the  three  groups,  about  800  -  900.  The  spatial  inhomogeneity  in 
attenuation  quality  factor  Q  near  station  WMQ  shows  good  correlation  with  the  tectonic 
structure  in  the  region. 


INTRODUCTION 
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The  Hen  Shan  Range  extends  east-west  for  at  least  2500  km.  This  large  tectonic 
belt  in  Central  Asia  lies  just  behind  the  collision  of  India  into  the  Asia  plate  and  forms  the 
most  active  intracontinental  mountain  range  in  the  world.  This  region  is  very  complex 
both  from  geology  and  topography,  and  there  are  many  active  thrust  and  strike-slip  faults 
(Tappomier  and  Molnar,  1979).  As  a  result,  the  seismic  properties  in  the  region  will  have 
large  spatial  variation. 

From  the  perspective  of  western  scientists,  the  region  is  remote  and  has  not  been 
thoroughly  studied.  But  it  has  become  interesting  to  seismologists  because  of  tectonic 
and  political  reasons.  The  development  of  methodology  makes  it  possible  to  study  the 
Hen  Shan  Range  seismologically  even  from  single  station. 

The  seismographic  station  WMQ  is  located  within  this  complex  region,  at  the 
north  boundary  of  Tien  Shan.  The  available  data  from  WMQ  in  the  late  1980’s  to  the  mid 
1990’s  consist  of  three  component  continuously  recorded  broadband  and  triggered  short 
period  digital  records,  which  are  well  suited  for  recording  regional  seismic  waves.  We 
selected  the  period  1988-1989  because  we  have  good  bulletin  information  in  these  years, 
from  Chinese  seismologists,  and  hence  have  epicenter  parameters  for  the  sources  of  our 
digital  signals. 

Patton  etal.  (1985)  studied  the  regional  structure  in  the  path  from  East  Kazakhstan 
to  the  Dzhungaria  basin  which  is  north  west  of  WMQ  and  Tien  Shan.  The  structure  in 
Tien  Shan  and  the  Tarim  basin,  to  the  south  of  Hen  Shan  are  likely  to  show  significant  dif¬ 
ferences  because  of  the  effects  of  the  collision  between  Indian  and  Asian  plates.  Gao  and 
Richards  (1995)  presented  evidence  that  the  depth  of  Moho  in  Xinjiang  is  quite  deep  (46 
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kilometers),  and  that  there  is  a  lateral  velocity  inhomogeneity  near  the  WMQ  station  based 
on  the  analysis  of  a  polarization  anomaly  from  P  waves  from  teleseismic  sources. 

In  this  paper,  we  shall  focus  on  the  spatial  attenuation  variation  of  shear  waves  in 
the  region  near  WMQ.  Our  overall  goal  is  to  achieve  a  better  understanding  of  attenua¬ 
tion,  in  a  detailed  study  of  the  propagation  of  seismic  waves  to  assist  in  the  interpretation 
of  seismic  energy  release  in  this  region. 

ATTENUATION  ESTIMATION 

Seismic  attenuation  represents  the  loss  of  seismic  energy  as  a  wave  propagates, 
over  and  above  what  is  expected  from  geometrical  spreading.  Attenuation  is  usually  fre¬ 
quency  dependent  especially  for  bands  that  include  high  frequency  signals. 

From  discussions  by  Atkinson  (1992)  and  Shi  etal.  (1996),  the  geometrical  spread¬ 
ing  of  Lg  waves  at  short  distances  is  itself  complex  and  does  not  follow  theory  such  as  the 
[A^^sinA^^]'^  dependence  predicted  for  Airy  phases.  The  direct  inversion  of  attenuation 
quality  factors  from  regional  earthquakes  is  therefore  difficult  and  probably  requires  the 
use  of  laterally  varying  structures. 

However,  the  empirical  Green’s  function  method  makes  it  possible  for  us  to  first 
obtain  the  source  information  and  then  to  estimate  the  path-dependent  attenuation  factor  Q 
from  fitting  the  Lg  displacement  spectrum  over  a  wide  frequency  range  by  using  the 
known  source  information.  This  process  does  not  require  the  geometrical  spreading  factor 
which  is  so  uncertain  for  regional  earthquakes. 
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We  assume  the  frequency  dependent  attenuation  spectrum  along  the  path  from 
event  i  to  station  j  as. 


with  A  being  the  epicentral  distance,  v  the  S-wave  velocity,  and  t  the  travel  time.  The 
attenuation  quality  factor  Q  is  usually  considered  to  be  frequency  dependent  {e.g.  Singh 
and  Herrmann,  1983;  Campillo,  1987)  and  we  assumed  Q  as, 

G  =  , 

with  Qo  =  Ihe  quality  factor  value  at  1  Hz,  and  q  =  the  exponent  parameter  of  the  quality 
factor  which  is  constant  over  the  frequency  band  of  interest  (e.g.  Mitchell,  1980). 

In  order  to  obtain  an  accurate  estimate  of  the  two  unknown  parameters  in  the  atten¬ 
uation  quality  factor,  we  must  first  eliminate  the  trade-off  problems  of  source  comer  fre¬ 
quency  and  Qo,  which  can  give  the  same  shape  of  the  signal  spectram  by  increasing  one 
parameter  and  decreasing  the  other.  However,  one  of  these  parameters  is  determined  by 
the  source  of  the  event,  and  the  other  parameter  is  determined  by  the  propagation  path 
from  event  to  stations.  Based  on  this  difference,  we  first  use  the  empirical  Green’s  func¬ 
tion  deconvolution  process  to  obtain  the  source  information,  especially  the  comer  fre¬ 
quency  of  the  event.  Then  we  introduce  the  corner  frequency  estimated  from  empirical 
Green’s  function  deconvolution  into  the  process  of  attenuation  analysis.  As  shown  by  Shi 
et  al.  (1996),  this  two  step  analysis  makes  it  possible  to  measure  the  two  spectral  parame¬ 
ters  accurately,  and  give  good  information  for  the  path  dependent  spatial  variation  of  the 
attenuation  quality  factor. 
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In  the  empirical  Green’s  function  deconvolution  process,  we  require  an  event  pair 
with  magnitude  different  enough  but  having  the  same  wave  propagation  path  and  same 
source  mechanism,  so  that  the  smaller  event  can  be  treated  as  the  empirical  Green’s  func¬ 
tion  of  the  larger  event  and  the  deconvolution  of  the  same  phases  of  the  two  events  will 
only  contain  source  information.  From  the  deconvolution,  which  is  a  spectral  ratio  in  the 
frequency  domain,  we  can  estimate  the  comer  frequency  of  the  larger  event  in  the  pair 
using  an  optimizing  method.  In  order  to  simplify  the  problem,  we  assume  the  source  has 
only  one  comer  frequency,  which  gives  an  upper  limit  for  the  size  of  the  larger  event  in 
each  pair.  Also  the  signal  to  noise  ratio  of  the  smaller  event  needs  to  be  good  enough  in 
the  frequency  band  for  analysis,  and  this  gives  a  limit  on  the  epicentral  distance  of  the 
events  and  a  lower  limit  on  the  size  of  the  smaller  event  in  the  pair. 

The  observed  spectmm  of  event  i  at  station  j  can  be  modelled  as, 

with  source  spectrum 

Here  Qqj  the  low  frequency  spectral  asymptote  and  is  proportional  to  the  scalar  seismic 
moment.  is  the  radiation  pattern  parameter,  and  is  the  source  corner  frequency 
based  on  Aki’s  (1967)  model.  Parameter  1(f)  represents  the  instrument  response,  F(f) 
the  site  response,  R  the  geometrical  spreading,  and  Bjj  the  event  to  station  path-averaged 
attenuation  term  in  frequency  domain. 
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The  factors,  /(/).  F(f),  R  and  By  cancel  out  in  the  spectral  ratio  we  are  considering, 

so  that 


r{f) 


1 

1  + v//,,)^ 


'f^Y 


Cq  is  the  ratio  of  Qq!  ^2- 

We  use  an  optimization  method  in  which  comer  frequencies and  j^2  varied 
on  a  grid  to  minimize  the  root  mean  square  of  -  rif)  and  to  find  the  best  fit^i, 

where  rQg^  is  the  spectral  ratio  from  an  observed  data  pair. 

The  second  step  is  to  estimate  the  two  parameters  Qq  and  t|  from  the  S  and  Lg  dis¬ 
placement  spectmm  of  the  signal.  We  define  a  prediction  root  mean  square  error. 


(  -1  ^ 

• 

'1 

ln[Ayif)  Ij  C/)J- 

In 

V 

Uq„/"J 

-c 

where  II II  denotes  the  root  mean  square  error  and  C  =  logg(Q.Qj-RiY-Fj)  which  is  considered 
as  frequency  independent.  Performing  a  grid  search  over  Qq  and  t|,  with  the  above  C 
being  the  one  that  minimizes  E,  we  choose  the  best  fit  of  Qq  and  T|  for  each  path,  given 
comer  frequency,  of  the  event.  When  we  assume  Fj  as  frequency  independent,  we  are 
ignoring  the  site  response  which  can  be  a  significant  factor  (e.g.,  Chun  et  al,  1987).  As  a 
result,  the  variation  in  the  attenuation  quality  factor  Q  that  we  have  estimated  includes  the 
site  effect  near  the  source  and  receiver  region. 


DATA  PROCESSING 
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The  seismograms  used  in  this  study  were  recorded  by  CDSN  (Chinese  Digital 
Seismographic  Network)  station,  WMQ  (Urumqi).  There  are  two  instrument  responses, 
one  flat  to  velocity  from  0.05  to  4  Hz  and  recorded  at  20  sample/s,  another  flat  to  velocity 
from  1  to  8  Hz  and  recorded  at  40  sample/s. 

We  selected  41  event  pairs  from  1988  to  1989  based  on  the  requirements  of  the 
empirical  Green’s  function  method.  The  station  WMQ  is  the  nearest  digital  station  to  the 
Lop  Nor  nuclear  test  site  in  Xinjiang,  China.  The  events  we  selected  (Table  1)  are  all  nat¬ 
ural  earthquakes  that  happened  quite  near  to  the  test  site,  with  paths  to  WMQ  crossing 
Tien  Shan  Mountains.  The  study  of  these  events  is  very  useful  to  analyzing  the  structure 
in  this  region.  Each  pair  of  events  selected  have  similar  foci  location,  similar  source 
mechanism,  and  the  waves  from  both  events  traveled  a  similar  path  and  are  recorded  by 
the  WMQ  station  with  similar  wave  forms,  but  with  different  magnitude  (Figure  1).  Some 
pairs  have  the  same  larger  events  and  there  are  only  28  independent  larger  events  in  the 
total  set  of  pairs.  The  pair  events  generally  have  focal  depth  within  0-20  km  range,  mag¬ 
nitude  from  2.7  -  5.1  magnitude  for  larger  events,  and  2.5  -  3.5  ;n^  for  smaller  events 
(Table  1).  The  spatial  distribution  of  these  events  are  mainly  south  and  southwest  of  the 
WMQ  station,  ranging  from  90  -  260  km  epicentral  distance. 

We  computed  the  displacement  amplitude  spectrum  of  vertical  component  P 
waves  using  about  2.0  seconds  window  after  the  first  P  anival  (mainly  Pg),  and  vertical 
component  S  and  Lg  waves  within  a  group  velocity  window  from  3.0  to  3.6  km/s.  The 
signal  to  noise  ratio  for  both  events  in  each  pair  is  well  above  2.  The  instrument  response 
was  removed  before  calculating  the  displacement  spectrum. 
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From  simplicity  in  analysis,  we  assumed  the  source  for  both  events  in  the  pairs  con 
forms  to  Aki’s  model  with  one  comer  frequency  for  the  source  spectram,  and  we 
ignored  the  directivity  effect  of  the  source  on  comer  frequency  measurement.  Discrepan¬ 
cies  among  estimates  for  large  local  events  may  come  from  violations  of  the  point  source 
approximation  and  possible  directivity  of  the  source,  and  one  station  estimation  of/c  can 
have  large  bias.  However,  for  earthquakes  with  magnitude  less  than  5,  the  duration  of  the 
source  is  only  about  1-2  seconds,  therefore  the  source  may  not  be  so  complex  as  having 
two  comer  frequencies  in  its  spectmm,  and  the  directivity  is  not  so  significant.  From  our 
previous  work  in  Northeastern  United  States  earthquakes  {Shi  et  al.  1998),  the  uncertainty 
caused  by  these  simplifications  is  not  larger  than  2  Hz  in  comer  frequency  estimation. 

Figure  2  is  an  example  of  our  procedure  for  estimating  the  comer  frequency  of  the 
larger  event  in  the  pair  by  searching  the  optimal  corner  frequency  and  f(.2  with 
source  model  and  frequency  dependent  Q  simulation.  Usually  we  only  estimated  comer 
frequency  for  the  larger  event  in  the  pair  and  fixed  the  corner  frequency  of  the  smaller 
event /^2  some  value.  The  uncertainty  in/^2  effect  the  estimation  of  ff-i  if  the 

size  difference  of  the  two  events  is  not  less  than  about  one  magnitude  unit. 

By  simulating  the  spectral  ratio  of  S  and  Lg  waves  within  a  certain  group  velocity 
window  or  first  several  seconds  of  P  waves  in  each  pair,  we  estimated  the  comer  frequen¬ 
cies  of  28  larger  events  in  the  pairs  and  plotted  the  estimated  comer  frequencies  with  the 
body  wave  magnitude  m/,  of  these  28  events  (Figure  3).  As  shown  in  the  plot  and  as 
expected,  the  comer  frequencies  for  these  events  decrease  as  the  size  of  the  earthquakes 
increase.  We  need  to  note  that  the  earthquakes  we  studied  here  all  have  size  larger  than 
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magnitude  3  except  one  event.  The  comer  frequency  relation  with  event  size  here  does 
not  conflict  with  our  previous  result  (Shi  et  aL,  1998)  from  earthquakes  in  Northeastern 
U.S.  that,  for  intraplate  earthquakes,  the  comer  frequency  of  the  event  does  not  depend  on 
the  event  size  when  the  earthquake  is  less  than  about  magnitude  3,  and  the  comer  fre¬ 
quency  decreases  with  increasing  event  size  for  larger  events. 

From  our  experience,  the  trade-off  between  comer  frequency  and  attenuation  qual¬ 
ity  factor  Q  is  the  main  problem  in  measuring  the  two  parameters  in  quality  factor,  Qq  and 
r\.  By  using  empirical  Green’s  function  deconvolution  method  in  comer  frequency  esti¬ 
mate,  we  measured  the  source  comer  frequency  from  the  source  spectmm  first  before 
introducing  the  comer  frequency  measurement  into  the  attenuation  estimation  process. 
The  trade-offs  are  then  sufficiently  eliminated.  Figure  4  shows  two  examples  in  the  simu¬ 
lation  of  the  source  and  attenuation  spectmm  model  with  the  S  and  Lg  displacement  spec¬ 
tmm  to  find  the  best  fit  two  path-averaged  attenuation  parameters  Qq  and  t|.  From  each 
spectmm,  we  estimate  an  event  -  station  path-averaged  Qq  and  q  and  obtain  the  spatial 
variation  of  attenuation  in  the  region  we  analyzed. 

Since  only  28  events  have  their  source  corner  frequency  well  estimated  from 
empirical  Green’s  function  analysis,  we  obtain  the  event-station  values  of  Qq  and  q  for  28 
paths  based  on  the  displacement  spectmm  of  S  and  Lg  waves  within  the  group  velocity 
range  from  3.6  to  3.0  km/s.  The  path  dependent  values  of  these  two  parameters  are  plotted 
in  Figure  5a  and  5b.  The  degree  of  the  background  shade  represents  the  elevation  in  the 
region,  with  mountains  showing  whiter. 

It  is  very  obvious  that,  for  the  quality  factor  Qq,  the  values  from  the  paths  which 
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cross  the  Tien  Shan  mountains  southwestern  of  WMQ  are  lowest,  about  500  to  600.  For 
the  paths  along  the  Hen  Shan  mountains  to  the  west  of  WMQ,  the  Qq  values  are  relatively 
higher,  about  700.  For  the  paths  through  the  basin  to  the  southeast  of  WMQ,  Qq  values  are 
also  relatively  higher,  about  800  -  900.  There  is  only  one  path  that  crosses  through  the 
Tien  Shan  on  east  side  of  WMQ,  and  Qq  seems  lower  than  the  path  values  of  other  groups, 
about  500.  Generally,  in  each  region,  Qq  from  S  and  Lg  waves  on  different  paths  are  found 
to  be  consistent  for  regions  of  similar  topography. 

The  variation  of  rj  for  different  paths  seems  not  so  much  structure  related  (as  is  the 
case  with  Qq),  but  is  more  randomly  variable. 

DISCUSSION  AND  CONCLUSIONS 

We  have  analyzed  Lg  spectra  from  28  earthquakes  in  XinJiang,  China,  with  body- 
wave  magnitude  ranging  between  2.7  to  5.1.  The  empirical  Green’s  function  method 
was  used  to  determine  the  source  corner  frequencies,  and  the  spectral  fitting  of  Lg  dis¬ 
placement  spectra  was  used  to  determine  the  path  averaged  Qq  and  t|  (Lg  2  at  1  Hz  and  its 
power-law  frequency  dependent  parameter,  respectively).  This  method  efficiently  elimi¬ 
nated  the  trade-off  between  source  corner  frequency  and  attenuation  quality  factor  Q  in 
the  estimate. 

The  Tien  Shan  is  situated  north  of  the  Tarim  basin,  consisting  of  approximately 
east-west-trending  blocks  which  are  often  separated  by  wide  latitudinal  valleys  (Ni,  1978). 
Abundant  evidence  for  thrust  faulting  throughout  the  present  Tien  Shan  indicates  that  the 
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present  range  was  built  by  crustal  shortening  and  crustal  thickening.  It  is  one  of  the  most 
seismically  active  mountain  chains  in  the  world,  with  lateral  variations  related  to  Ceno- 
zoic  crustal  shortening  {Burov  et  al,  1990;  Li,  1990).  Hence,  the  geological  structure  near 
WMQ  is  quite  complex  as  shown  by  Feng  (1985),  and  Patton  et  al.  (1985).  We  will 
expect  the  spatial  variation  in  attenuation  corresponding  with  this  structure  will  also  be 
complex.  In  order  to  see  the  correlation  of  the  attenuation  we  obtained  with  the  local 
structure,  we  contoured  the  values  of  Qq  based  on  the  path  averaged  Qq  information  from 
S  and  Lg  waves  (Figure  6).  The  region  we  analyzed  can  be  considered  as  three  groups, 
one  group  along  the  high  elevated  Tien  Shan,  another  group  cross  the  Tien  Shan  suture 
zone  to  southwest  of  WMQ,  and  one  within  the  relatively  lower  elevated  basin  region  to 
the  southeast  of  WMQ. 

The  spatial  inhomogeneity  in  attenuation  near  Urumqi  is  clearly  correlated  with 
the  tectonic  structure  and  path  direction  in  the  region.  For  the  waves  propagating  along 
the  Tien  Shan  mountains  west  of  WMQ  and  south  into  the  basin  region  which  do  not  cross 
the  thrust  fault  zone  shown  by  Tapponnier  (1979)  and  Feng  (1985),  the  quality  factor  Q  is 
relatively  higher,  and  there  is  relatively  lower  attenuation.  While  for  the  seismic  waves 
propagating  across  the  Hen  Shan  mountains  to  the  southwest  which  cross  the  thrust  fault 
zone,  we  found  the  quality  factor  Q  is  relatively  lower. 

The  attenuation  we  measured  is  an  apparent  attenuation  which  includes  both 
intrinsic  and  scattering  phenomena  within  the  upper  crust  in  the  area.  The  quality  factor  Q 
at  different  frequencies  is  actually  the  total  effect  from  intrinsic  and  scattering  attenuation 
based  on  the  scale  of  the  structure  at  various  frequencies. 
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HGURE  CAPTIONS 

Figure  1.  An  example  of  event  pairs  in  January  9,  1989  (mi,  =  3.8  and  2.8) 
recorded  by  same  station  WMQ.  The  amplitude  of  the  smaller  event  {OT:  01/09/1989 
15:22:01)  is  almost  10  times  smaller  than  the  larger  event  (07:  01/09/1989 15:10:23)  but 
with  similar  wave  forms. 

Figure  2.  Examples  showing  the  spectral  ratio  (thin  lines)  of  two  event  pairs  from 
Lg  waves  within  a  group  velocity  window  from  3.0  to  3.66  km/s  and  P  waves  of  first  2.0 
seconds  after  first  P  arrival.  The  thick  smooth  curve  is  the  synthetic  source  function  ratio 
in  the  frequency  domain  with  estimated  comer  frequency  information  from  optimization 
method.  The  solid  arrow  shows  the  position  of  the  comer  frequencies  of  the  larger  event 
(fci)  and  the  smaller  event  (^^2)-  If/c2  is  not  estimated,  we  fixed /c2  =  40  Hz  and  use  an 
empty  arrow  to  show  it. 

Figure  3.  Two  examples  of  the  displacement  spectra  from  which  we  estimated  the 
attenuation  information,  Qq  and  Tj,  with  calculated  corner  frequency  information.  The 
arrows  represent  the  position  of  the  corner  frequency  of  the  event. 

Figure  4.  Plot  of  corner  frequency  with  magnitude  of  the  larger  events  of  the 
pairs.  The  squares  represent  corner  frequencies  from  the  P  wave  spectral  ratio.  The  trian¬ 
gles  represent  comer  frequencies  estimated  from  the  Lg  wave  spectral  ratio. 

Figure  5a.  Qq  values  for  all  the  28  event-station  paths  we  analyzed  from  Lg  waves. 

Figure  5b.  rj  values  for  the  event-station  paths. 

Figure  6.  Attenuation  quality  factor  Qq  contours  with  the  background  topography 
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in  the  region.  The  density  of  the  shade  shows  the  elevation  in  the  region,  with  whiter 
shade  being  higher  elevated. 
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Table  The  event  pairs  selected  from  earthquakes  in  XinJiang. 
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18 
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04 
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03 
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a 
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02 
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11 
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- 
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13 
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88.433 

- 

2.5 

3 
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04 

06 

21 

15 

23.8 

42.567 

85.833 

14. 

2.7 

9.0 

8.0 

1988 

01 

05 

14 

14 

37.3 

42.617 

85.850 

. 

2.6 

4 
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07 

04 

00 

54 

34.2 

43.683 

87.%7 

28. 

4.6 

3.25 

2.5 

1988 

07 

04 

01 

01 

24.2 

43.683 

87.967 

26. 

3.4 

5 

1988 

08 

03 

23 

51 

47. 

43.633 

85.350 

9.0 

3.2 

_ 

11.5 

1988 

12 

30 

20 

41 

38.2 

43.633 

85.333 

12. 

2.5 

6 

1988 

08 

09 

02 

51 

23.9 

43.783 

86.533 

15. 

3.3 

. 

6.0 

1989 

11 

07 

12 

51 

31.1 

43.700 

86.617 

15. 

2.5 

7 

1988 

09 

03 

03 

50 

8.9 

43.333 

85.350 

18. 

3.6 

3.25 

3.0 

1989 

04 

26 

05 

47 

54.9 

43.017 

85.367 

18. 

2.6 

8 

1988 

09 

18 

06 

39 

11.0 

41.983 

88.683 

11. 

3.6 

6.75 

6.0 

a 

1989 

11 

28 

08 

37 

25.4 

41.983 

88.800 

16. 

2.8 

b 

1989 

03 

01 

16 

57 

9.9 

42.050 

88.550 

11. 

2.7 

9 

1988 

09 

24 

11 

12 

38.3 

43.367 

88.783 

13. 

3.7 

6.5 

1989 

04 

19 

04 

46 

12.6 

43.717 

88.583 

13. 

2.5 

10 
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10 

05 

17 

10 

22.9 

42.417 

85.850 

23. 

3.4 
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a 
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06 

26 

20 

43 
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14 
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06 

06 

01 
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12 

21 

16 

22 
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10 

10 
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2. 
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18 
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11 
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15 
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11 

15 

16 

56 

45.3 

42.217 

89.333 
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5.1 

2.33 

1.67 

a 
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11 

15 

16 

54 
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42.200 

89.400 
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2.5 
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07 

23 

15 
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15 
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12 

07 

05 

58 
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Note, 


The  first  event  in  each  pair  is  the  larger  event,  and  the  smaller  ones  are  the  possible 
pair  event  to  the  larger  event.  Some  event  may  have  more  than  one  possible  pair,  and  the 
comer  frequencies  we  estimated  are  the  average  of  values  derived  by  deconvolution  of  all 
possible  pairs. 
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Part  3.  Description  of  an  unusual  swarm  of  seismic  activity  in  Western  China 

The  Lop  Nor  test  site,  like  the  Nevada  test  site  of  the  United  States,  lies  in  an  area  of  active 
tectonics,  and  thus  will  be  associated  with  numerous  seismic  events  that  may  be  detected  and  are 
likely  to  require  analysis  in  a  CTBT  monitoring  program.  This  section  of  the  present  report 
describes  an  unusual  swarm  of  seismic  activity  that  is  a  subset  of  the  Lop  Nor  seismicity  described 
in  a  previous  report  under  our  grant  from  AFOSR.  The  swarm  consisted  of  more  than  400  small 
seismic  events  in  a  period  of  about  100  days,  beginning  on  December  27, 1975,  and  centered  in 
Western  China  at  (44.25°N,  85.65°E).  It  is  not  clear  whether  the  swarm  was  a  natural 
phenomenon,  or  whether  it  was  due  to  human  activity  (such  as  mine  blasting). 

To  give  the  context.  Figures  1  through  5  of  this  part  of  the  final  report  show  the  location  of 
2757  seismic  events  reported  by  the  Chinese  State  Seismological  Bureau,  for  the  years  1970  and 
1972  through  1989,  for  the  region  between  latitudes  37°N  and  45°N,  and  between  longitudes  85°E 
and  93®E.  Figure  1  shows  the  location  of  all  the  2757  events,  which  have  magnitude  2  and  above 
(the  full  set  of  events),  together  with  a  rectangle  close  to  (41.5°N,  88.5°E)  indicating  the  location  of 
the  Lop  Nor  nuclear  test  site.  Figure  2  shows  the  same  data  set  for  the  884  events  of  magnitude  3 
and  larger;  Figure  3  is  for  the  164  events  of  magnitude  4  and  larger;  Figure  4  is  for  the  31  events  of 
magnitude  5  and  larger;  and  Figure  5  is  for  the  4  events  of  magnitude  6  and  larger. 

When  the  events  for  1975  and  1976  are  examined  for  their  distribution  in  time,  in  space, 
and  in  magnitude,  as  shown  in  Figures  6  and  7,  it  is  apparent  that  several  hundred  events  are 
strongly  clustered  in  time  and  space,  beginning  on  December  27, 1975,  and  ending  about  one 
hundred  days  later.  These  434  events  are  all  within  20  km  of  the  point  (44.25°N,  85.65°E),  and 
about  300  km  Northwest  of  the  Lop  Nor  test  site. 

Figure  8  shows  on  an  expanded  scale  the  location  of  the  434  events  occurring  in  the  swarm 
from  Dec  27, 1975,  to  April  4, 1976.  It  is  apparent  that  the  State  Seismological  Bureau  reports 
event  locations  to  the  nearest  1/60  of  a  degree  (i.e.,  to  the  nearest  minute),  so  the  locations  are 
discretized  and  many  of  the  points  shown  represent  multiple  events.  The  Times  Atlas  shows  an 
East-West  road  running  at  about  latitude  44.2°N  through  the  points  in  this  Figure,  and  to  the  town 
of  Manass  off  the  Figure  to  the  East  at  (44.20°N,  86.14°E). 

Figure  9  shows  the  number  of  events  reported  per  day.  There  is  no  particular  pattern,  as 
for  example  might  expected  for  mine  activity  (for  which  a  weekly  periodicity  is  common,  some 
mines  choosing  to  blast  mainly  on  weekends,  and  other  mines  avoiding  weekends).  Figure  10 
shows  the  events  occurring  per  half-hour  of  local  time-of-day,  for  the  100-day  duration  of  the 
swarm,  and  again  there  appears  to  be  no  simple  pattern  such  as  a  dominance  of  events  during 
daylight  working  hours  (very  common  for  mining  activity  and  construction  blasting). 

It  is  common  for  an  earthquake  to  be  followed  by  a  sequence  of  smaller  events.  In  such 
cases  the  sequence  begins  with  a  large  event  (the  "main  shock"),  followed  by  aftershocks  of 
general  decreasing  magnitude  and  frequency.  Figure  1 1  shows  how  the  434  events  are  distributed 
in  magnitude  and  time.  Such  a  pattern  is  absent  in  this  figure.  The  sequence  simply  ends  in  early 


April,  as  shown  also  in  Fig.  7.  For  a  swarm  of  earthquakes  associated  with  volcanic  activity,  there 
is  typically  no  dominant  event  that  starts  the  sequence — which  tends  to  have  a  gradual  beginning 
and  ending. 

Figure  12  shows  how  the  434  events  are  distributed  in  magnitude.  In  this  plot  of  N(M) 
against  M,  N(M)  is  the  cumulative  number  of  events  in  the  cluster  with  magnitude  M  or  greater.  A 
line  given  by  logN(M)  =  5.4- 1.4  M  is  a  good  fit  to  the  data.  A  slope  closer  to  the  value 
-1  would  be  more  typical  of  a  natural  aftershock  sequence. 

The  swarm  is  unusual  in  that  it  does  not  display  the  characteristics  of  either  a  natural 
earthquake  sequence,  or  human  activity  such  as  blasting  associated  with  a  construction  project  or 
mining.  It  is  presumably  one  or  other  of  these  alternatives.  Possibly,  the  sequence  continued  but  a 
decision  was  made  not  include  the  events  following  April  4, 1976  in  the  State  Seismological 
Bureau  catalog.  Figures  6  and  7  indicate  that  only  for  the  duration  of  the  swarm,  were  events 
reported  that  had  magnitude  lower  than  2.5.  Since  the  daily  rate  often  reached  around  10  or  20 
events,  and  occasionally  even  more  (see  Fig.  9),  mine-  or  constraction-blasting  seems  unlikely.  It 
cannot,  however,  be  ruled  out.  If  a  swarm  this  large  raised  questions  in  the  context  of  CTBT 
monitoring,  the  nature  of  the  swarm  could  perhaps  best  be  resolved  by  non-seismic  methods. 
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Figure  1.  The  location  of  all  2757  seismic  events  for  the  19-year  period  1970,  and  1971  -  1989, 
reported  by  the  Chinese  State  Seismological  Bureau  for  the  rectangle  from  37°^5°N,  85°-93°E. 
The  location  of  the  Lop  Nor  test  site  is  shown  as  the  rectangle  41.3°-41.8°N,  88.3°-88.8°E. 


Figure  2.  Same  as  Fig.  1,  but  showing  the  884  events  with  magnitude  3  and  larger.  Note  a 
concentration  of  events  near  (44.25°N,  85.65°E).  These  events  clustered  in  time  as  well  as  space 


Figure  3.  Same  as  Fig.l,  but  showing  the  164  events  with  magnitude  4  and  larger. 


Figure  4.  Same  as  Fig.l,  but  showing  the  31  events  with  magnitude  5  and  larger. 


Figure  6.  Distribution  of  seismicity  in  the  Lop  Nor  region,  with  respecUo  magnitude,  time,  and  space,  for  the  year  1975 
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Figure  7.  Distribution  of  seismicity  in  the  Lop  Nor  region,  with  respect  to  magnitude,  time,  and  space,  for  the  year  1976 


Figure  8.  Map  of  locations  near  (44.25°N,  85.65°E)  of  434  clustered  events  in  W.  China 
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Figure  10.  The  number  of  events  per  half-hour  is  shown  as  a  function  of  local  time-of  day.  No 
pattern  is  apparent,  such  as  occurrence  during  daylight  hours  or  the  usual  working  hours. 


cluster,  W.  China;  100  days;  log  N(M)  =  5.4  -  1.4  M 


M;  near  (44.25°N,  85.65°E);  27  Dec  '75  to  4  Apr  '76 


Figure  12.  The  cumulative  magnitude,  N  =  N(M)  =  number  of  events  with  magnitude  >  M,  is 
shown  to  fit  a  straight  line  of  the  form  log  N(M)  =  5.4  -  1.4  M.  Such  a  line  is  common  for 
earthquake  activity,  but  not  for  a  cluster  of  short  duration  that  does  not  begin  with  a  main  shock  of 
larger  magnitude. 
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ABSTRACT 

The  use  of  seismic  arrays  to  determine  slowness  and  travel  time  greatly  enhances  the 
identification  of  later  arrival  branches  of  the  P-wave  travel  time  curve  between  12®  and  30®.  The 
cross-array  in  Kurchatov,  Kazakstan  is  a  medium  aperture  array  with  a  total  aperture  length  of 
22.5  km.  It  has  provided  high-quality  digital  records  and  is  well  located  for  a  detailed 
investigation  of  travel  times  and  slownesses  of  P-wave  branches  at  distances  around  20® 
underneath  most  parts  of  central  Asia.  In  this  paper,  over  100  earthquakes  and  an  underground 
nuclear  explosion  which  occurred  in  the  Tibet  Plateau  and  western  China  during  1995-1996  were 
used  to  study  the  upper  mantle  velocity  structure.  Events  which  have  occurred  in  the  distance 
range  10®  -  32®  and  with  event-to-station  azimuth  range  of  312°  -  339°  from  the  Kurchatov  array 
are  selected  for  analysis.  A  simple  weighted,  delay-and-sum  beam-forming  method  was  used  to 
determine  the  travel  times,  slownesses  and  azimuths  of  the  first  arrival  P  as  well  as  later  arrival  P 
phases.  The  travel  time  and  slowness  data  are  corrected  for  focal  depth  and  local  structure,  and  are 
used  to  invert  for  the  upper  mantle  P-wave  velocity  structure.  The  results  of  the  inversion  for  one 
dimensional  velocity  structure  indicate  the  existence  of  several  discontinuities  including  known 
discontinuities  at  depth  of  417  km  and  690  km.  A  preliminary  result  is  that  a  sharp  increase  in  P 
velocity  between  about  550  to  600  km  depth  is  indicated  by  the  later-arrival  P  phases.  The  data 
provide  resolution  of  velocity  stmcture  down  to  about  780  km  depth. 
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Introduction 

With  the  advent  of  digital  seismometry  and  recent  advances  in  microprocessor  and 
communication  technologies,  there  is  an  upsurge  of  interest  in  seismic  arrays.  Approximately  22 
to  30  seismic  arrays*  are  being  considered  for  deployment  in  the  next  few  years  as  primary 
network  stations  of  the  International  Monitoring  System  (IMS)  for  monitoring  the  Comprehensive 
Test  Ban  Treaty.  Availability  of  the  seismic  data  from  these  IMS  seismic  arrays  deployed  at 
various  places  worldwide  will  certainly  stimulate  the  use  of  seismic  array  data  for  various 
seismological  investigations.  The  cross-array  in  Kurchatov,  Kazakstan  is  a  such  regional  medium 
aperture  array.  Figure  1  shows  the  cross-array  location  and  its  configuration.  The  large  territory 
of  the  central  Asian  Republic  of  Kazakstan,  located  in  the  stable  interior  of  Asia,  is  uniquely  suited 
for  sensitive  monitoring  of  global  seismicity,  and  nuclear  explosions.  The  great  distance  from 
major  oceans,  low  levels  of  cultural  and  seismic  activity,  and  the  excellent  transmission 
characteristics  of  the  exposed  rock  of  the  Kazakh  platform  provide  locations  for  low-noise,  low- 
attenuation  recording  sites.  In  October  of  1994,  a  centralized  digital  recording  system  was  installed 
to  record  the  short  period  seismometers  of  an  existing  21 -element  cross-array  at  Kurchatov  city  on 
the  former  East  Kazakstan  nuclear  test  site.  The  cross-array  has  two  orthogonal  linear  arms,  each 
consisting  of  1 1  boreholes  (average  depth  «  45  m)  containing  short-period,  vertical-component 
seismometers  with  station  spacing  of  2.25  km  (Kim  et  ai,  1996).  The  total  aperture  length  is 
about  22.5  km.  The  instrument  response  of  the  borehole  seismograph  is  nearly  flat  to  input 
ground  velocity  in  the  frequency  band  0.5  -  5  Hz  and  signals  are  recorded  continuously  with  a 
sampling  rate  of  40  samples/sec.  The  cross-array  has  provided  high-quality  digital  records  since 
its  operation  and  is  well  located  for  a  detailed  investigation  of  travel  times  and  slownesses  of  P- 
wave  branches  at  distances  around  20°  underneath  most  of  the  central  Asia. 

The  most  direct  information  on  the  P-wave  velocity  distribution  in  the  upper  800  km  of  the 
mantle  comes  from  observations  of  refracted  P-wave  arrivals  at  epicentral  distances  of  less  than 
30°.  Two  general  properties  of  the  velocity  structure  of  the  upper  mantle — a  zone  of  low  or 
negative  velocity  gradient  above  250  km  (the  low-velocity  zone),  and  regions  of  high  velocity 
gradient  between  300  and  700  km  (the  transition  zone) — are  now  accepted.  However,  regional 
models  differ  significantly  in  detail.  The  differences  may  be  caused  by  variations  in  structure,  but 
difficulties  in  the  interpretation  of  the  observed  data  are  also  involved. 

In  this  paper,  beam-forming  was  used  in  a  relatively  small  time  window  (usually  1 .2  to  2 
seconds)  to  determine  the  slownesses  and  travel  times  not  only  for  the  first  arrival  but  also  for  the 
later  arrivals.  Based  on  these  measurements,  we  derive  an  initial  upper  mantle  P  velocity  structure 

*  ANNEX  1  to  the  Protocol,  Comprehensive  Nuclear  Test-Ban  Treaty,  September  10,  1996,  listed  22  seismic 
arrays,  and  eight  additional  stations  which  will  each  be  deployed  as  a  three-component  station  to  be  upgraded  to  an 
array  at  a  later  date. 
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in  western  China  using  the  Herglotz-Wiechert  inversion  technique,  which  uses  distance  as  a 
function  of  dTIdA  .  We  then  modify  the  initial  model  by  comparing  the  predicted  travel  time  and 
the  observations. 


Data 

The  region  studied  is  one  of  the  most  complex  areas  of  present  day  tectonic  activity.  It 
appears  impossible  to  find  a  one-dimension  model  to  fit  all  data  from  a  wide  range  of  sources.  But 
in  a  narrow  azimuth  range  such  as  shown  in  the  Figure  2,  it  is  possible  to  find  a  one-dimension 
model  to  fit  the  observations  for  a  good  range  of  distances.  We  selected  approximately  100  events 
which  have  occurred  in  western  China  during  1995  and  the  first  four  months  in  1996.  The 
epicentral  distances  of  the  selected  events  extend  from  10°  to  beyond  30°  in  the  azimuth  range  312° 
to  339°  (from  the  event  to  the  cross-array)  to  cover  this  narrow  area.  The  events  are  listed  in  Table 
1  and  the  locations  of  epicenters  are  shown  in  Figure  2. 

A  seismic  record  section  of  the  waveforms  from  earthquakes  in  the  distance  range  10°  to 
20°  is  shown  in  Figure  3.  The  waveforms  recorded  at  the  center  of  the  array  (C21;  Figure  1)  are 
plotted  with  a  reduction  velocity  of  10  km/s  (=  1 1.12  s/deg).  The  record  at  a  distance  around  1 1.5° 
is  from  an  underground  nuclear  explosion  on  08/17/95  (00:59)  at  the  Lop  Nor  Chinese  test  site.  In 
Figure  3,  the  predicted  travel  time  curve  calculated  using  the  lASPEI  91  model  for  a  focal  depth  of 
15  km  is  plotted  and  various  P-wave  travel  time  branches  are  indicated  (Kennett  &  Engdahl, 
1991).  For  instance,  the  A-B  branch  represents  Pn  arrivals  and  the  C-D  branch  represents  first 
arrival  P-waves  whose  rays  have  turning  points  below  the  410  km  discontinuity,  whereas  the 
branch  G-H  is  a  portion  of  the  P  travel  time  branch  corresponding  to  rays  turning  below  the  660 
km  discontinuity.  A  line  joining  C  to  B  is  the  receding  branch  due  to  a  zone  of  rapid  velocity 
increase  near  the  410  km  discontinuity.  However,  the  data  do  not  show  a  clear  amplitude  change 
at  the  edge  of  the  receding  branch  at  C  in  Figure  3  (distance  range  of  14°-16°),  suggesting  a  less 
pronounced  velocity  change  above  the  410  km  discontinuity  in  the  study  area.  The  record  from 
each  event  is  plotted  without  adjusting  for  its  focal  depth.  Note  that  at  certain  distances  there  are 
multiple  arrivals  within  a  very  short  time  window,  as  well  as  possible  arrivals  of  depth  phases 
generated  near  the  source.  These  complex  arrivals  will  lead  to  slownesses  and  travel  times  that  are 
poorly  determined  for  parts  of  particular  travel  time  branches. 

A  synthetic  record  section  calculated  using  the  IAS  PEI  91  model  for  a  focal  depth  of  15  km 
is  shown  in  Figure  4  to  help  interpret  the  observed  record  section  shown  in  Figure  3.  The  source 
mechanism  from  an  event  on  12/18/95  (04:57,  M=5.7)  is  used  to  calculate  amplitude  and  travel 
time  using  the  WKBJ  seismogram  method  (Chapman,  1978).  The  synthetics  show  essentially 
similar  features  for  the  first  P  waves  when  compared  with  the  observed  record  section  shown  in 
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Figure  3.  However,  the  synthetics  show  less  complex  later  arrivals,  since  the  depth  phases  are  not 
included.  Also  simplifying  the  synthetics,  the  source  time  function  used  is  a  delta  function 
convolved  with  an  attenuation  operator  having  t*  =  0.1  s.  Nevertheless  the  synthetic  record 
section  is  a  useful  tool  to  interpret  various  later  arrivals  on  the  observed  records. 

Observations  of  Travel  Time  and  Slowness 

Analysis  of  Array  Data 

It  has  long  been  recognized  that  a  considerable  reduction  of  seismic  noise  is  possible  by 
beam-forming  array  data  (Mykkeltveit  et  al,  1990;  among  others).  We  employed  a  simple 
weighted,  delay-and-sum  (WDS)  beamforming  algorithm  to  process  the  cross-array  data  to 
measure  the  slowness,  p,  and  arrival  time,  T,  of  the  first  P  wave  as  well  as  several  later  arrival 
P  waves  at  the  center  of  the  array  (Capon  et  al.,  1967;  Johson  &  Dudgeon,  1993).  The  array 
beam,  z(t),  calculated  for  the  array  center  located  at  (xq,  yQ)  may  be  written  as 

Z(.t)  =  ^^w^^(t-6tm)  (1) 

-  ^o)  +  CVm  -  3^0)  +  V 

where  =  amplitude  weight  at  the  m-th  sensor  located  at  y^jj),  s^(t)  =  record  at  the 
m-th  sensor,  St^  =  time  delay,  </►  =  direction  of  the  arriving  plane  wave,  c  =  propagation 
velocity  of  the  plane  wave  (reciprocal  of  the  horizontal  slowness),  and  =  station  correction  that 
may  be  applied.  The  essence  of  the  method  is  that  each  record  is  shifted  in  time  with  respect  to  the 
reference  record  by  the  delay,  so  that  the  signal  arrival  time  becomes  common  to  all  stations;  and 
then  a  simple  algebraic  sum  of  all  the  records  is  computed  at  each  instance.  The  delays  compensate 
for  the  propagation  delays  and  the  waveforms  at  individual  sensors  add  constructively.  It  is  shown 
that  the  WDS  yields  nearly  as  much  signal  enhancement  as  other  more  elaborate  array  processing 
techniques  such  as  adaptive  array  processing  (e.g..  King  et  al.,  1973;  Filson,  1975;  Johnson  & 
Dudgeon,  1993).  In  equation  (1),  c,  (j)  and  are  assumed  to  be  known  quantities.  In  practice, 
wave  propagation  velocity,  c,  and  direction,  <j>,  are  found  by  calculating  all  beams  in  the  expected 
slowness  ranges,  say,  0.07  to  0.125  s/km  with  a  preset  increment  of  approximately  0.001  s/km  in 
both  X  and  y  directions,  and  then  finding  the  c  value  that  gives  the  clearest  beam  signal.  In  this 
study,  usually  over  1,000  beams  are  calculated  for  each  time  window  and  the  best  fit  is  found  that 
yields  maximum  beam  power,  which  is  computed  as  the  mean-squared  amplitude.  The  observed 
data  are  pre-filtered  with  a  third  order  Butterworth  filter  with  cutoff  frequencies  at  1  and  5  Hz. 

Usually  about  15  sec  of  signals  after  the  first  arrival  P-wave  are  analyzed.  Beam-forming 
is  performed  using  a  1  to  2  second  time  window  depending  upon  the  size  of  the  event,  since  large 
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events  tend  to  show  longer  period  pulses  due  to  source  duration.  We  find  strong  beam  power  in 
several  time  positions  which  indicate  the  arrivals  of  signals  when  we  move  the  beam-forming 
window  along  the  time  axis.  The  positions  of  high  beam  power  accurately  determine  the  arrivals 
of  signals. 

An  example  of  array  processing  of  the  data  is  shown  in  Figure  5.  The  signal  to  noise  ratio 
(SNR)  is  enhanced  by  the  beam-forming  technique,  and  the  identification  of  the  later  arrivals  is 
improved.  In  Figure  5a,  the  first  trace  is  the  reference  station  waveform,  and  the  second  trace 
shows  the  array  beam  calculated  by  using  the  slowness  and  azimuth  determined  from  the  small 
time  window  which  is  indicated  by  the  short  vertical  lines.  The  apparent  velocity,  azimuth,  and  the 
maximum  beam  power  determined  from  each  small  time  window,  are  plotted  in  Figure  5b.  When 
we  move  the  time  window  along  the  time  direction,  we  obtain  the  apparent  velocity,  azimuth  and 
beam  power  at  each  specific  time.  This  process  allows  us  to  identify  clearly  the  first  arrival  as  well 
as  depth  phases  and  other  later  arrivals.  It  is  clear  that  when  a  signal  arrives,  the  beam-power 
sharply  increases  from  the  background  level,  and  the  azimuth  and  slowness  become  stable.  We 
can  directly  obtain  the  slowness  and  travel  time  values  for  the  first  and  later  arrivals  by  carefully 
comparing  these  figures  (Figs.  5a  through  5d). 

First  Arrivals 

Usually  we  can  pick  up  the  first  arrivals  precisely.  Studying  them  give  us  some  important 
constraints  on  the  later  interpretation.  Figure  6  shows  the  travel  time  and  slowness  of  the  first 
arrivals.  Using  the  PDF  location  (Preliminary  Determination  of  Epicenters,  published  by  the  U.S. 
Geological  Survey)  for  the  source,  the  arrival  times  and  distances  are  corrected  to  surface  focus. 
The  measured  slowness  value  (1/c)  is  used  to  determine  the  corrections  for  each  arrival  using  the 
equations  given  by  Bullen  (1963,  p.  104).  For  the  sake  of  consistency  with  the  initial  model  used 
in  the  velocity  inversion  in  the  later  section,  the  lASPEI  91  velocity  model  has  been  used  in 
determining  these  corrections.  These  corrections  are  about  4.2  sec  and  0.13°  for  the  travel  time  and 
distance,  respectively,  for  a  P  arrival  with  a  slowness  of  0. 1  s/km  from  an  event  with  focal  depth 
of  20  km,  while  for  P  arrivals  with  the  same  slowness  from  an  event  with  depth  of  10  km,  these 
corrections  are  2.1  sec  and  0.06°  for  travel  time  and  distance,  respectively.  Hence,  even  a 
moderate  amount  of  error  in  the  event  focal  depth  will  affect  the  correct  interpretation  of  the  travel 
time  and  slowness  data.  The  focal  depth  is  usually  the  least  well  constrained  source  location 
parameter  in  the  PDE  listings.  For  most  of  the  events  with  unconstrained  focal  depth, 
conventionally  listed  in  the  PDE  catalog  as  33  km,  we  determined  the  focal  depth  using  slowness 
of  later  arrival  phases  and  we  also  adjusted  the  origin  times  of  these  events. 

When  the  travel  time  and  slowness  data  are  presented  separately,  the  relationship  between 
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travel  times  and  slownesses  is  not  readily  apparent,  and  resolution  of  finer  features  is  lost  in  the 
scatter  of  the  data.  In  Figure  6,  a  substantial  improvement  has  been  achieved  by  the  incorporation 
of  the  travel  time  and  slowness  data  in  a  single  presentation  following  the  method  of  Simpson  et 
al.  (1974).  At  each  travel  time  point  a  line  is  drawn  with  the  slope  corresponding  to  the  measured 
value  of  slowness  reduced  to  a  reference  slowness  of  0. 1  s/deg.  The  length  of  the  line  is  scaled  to 
the  weight  given  to  the  slowness  determination  and  corresponds  to  ±0.5®  for  a  well  determined 
slowness.  Part  of  the  scatter  in  the  data  shown  here  can  be  attributed  to  errors  in  the  PDE 
locations,  in  particular,  errors  in  focal  depth.  The  travel  time  and  slowness  data  are  also  subjected 
to  the  effects  of  local  structure  beneath  the  array  and  to  possible  systematic  bias  due  to  errors  in  the 
travel  time  curve  used  to  locate  the  events. 

A  careful  inspection  of  Figure  6  reveals  that  the  slope  lines,  which  should  lie  tangent  to  the 
trend  of  the  travel  times  at  each  point,  show  a  bias  away  from  the  average  curve.  This  bias  is  due 
to  a  shift  in  the  measured  values  of  slowness  caused  by  local  structure  near  the  array.  If  a 
knowledge  of  the  structure  causing  these  variations  is  available,  it  is  possible  to  compute 
corrections  to  the  measured  slowness  value  (Niazi,  1966;  Zengeni,  1970).  Niazi  has  shown  that 
the  effect  of  dipping  layers  beneath  an  array  is  to  ehange  the  measured  values  of  slowness  by  a 
factor  which  remains  essentially  constant  over  limited  ranges  of  azimuth  and  slowness.  Simpson 
et  al  (1974)  used  this  method  in  analyzing  the  WRA  (Warramunga,  Australia)  array  data.  We 
applied  this  method  to  the  Kurchatov  cross-array  with  events  in  a  narrow  azimuth  range.  At  each 
distance  where  a  slowness  determination  was  available,  the  slope  of  the  smoothed  curve  was 
compared  with  the  measured  value  of  slowness.  The  rms  residual  was  computed  and  the 
procedure  was  repeated  with  all  measured  first  arrival  slowness  values  multiplied  by  a  factor 
ranging  from  1.2  to  0.9.  A  minimum  in  the  rms  residual  occurs  with  a  correction  factor  of  1.06, 
indicating  that  the  observed  slowness  values  are  low  by  about  6%.  The  travel  times  with  slowness 
values  corrected  for  local  structure  are  shown  in  Figure  6,  which  shows  the  improvement  in  the 
trend  of  the  slope  lines  with  respect  to  the  travel  time  curve,  especially  between  20®  and  24°.  The 
abrupt  changes  in  slowness  in  the  first  arrival  data  of  Figure  6b  indicate  that  there  are  at  least  three 
major  first  arrival  branehes  with  cross-over  points  at  19®  and  25®. 

Later  Arrivals 

Based  on  the  previous  work,  we  can  pick  up  not  only  the  first  arrival  but  also  the  later 
arrivals  in  the  improved  waveform.  All  data  are  corrected  for  event  depth  and  local  structure. 
Figure  7  show  the  slowness  distribution,  and  the  size  of  the  circles  represent  the  quality  of  the 
slowness  measurement.  Even  though  the  scatter  in  the  data  cannot  give  us  clear  choices  of 
slowness  change,  we  can  still  recognize  that  there  exist  several  branches  in  order  to  fit  the 
distribution. 
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Velocity  Structure 


Herglotz-Wiechert  Inversion 

If  the  source  and  the  receiver  are  both  at  the  Earth's  surface,  the  angular  distance  A(p) 
between  the  source  and  the  point  at  which  the  ray  with  parameter  p  emerges  can  be  written  as  (see, 
e.g.,  Aki  and  Richards,  1980) 


A{p)  =  2  p 


dr 

r 


(2) 


where  p  is  the  ray  parameter,  p  =  rsini/  v(r),  rg  is  the  Earth's  radius,  r^  is  the  radial 
distance  to  the  deepest  point  of  the  ray  and  v  =  v(r)  is  the  velocity  as  a  function  of  radius.  The 
corresponding  travel  time  T(p )  is  given  by 


r^o 


T(p)  =  2 


\rlvy 


slirlvf-p^ 


dL 

r 


(3) 


The  Herglotz  and  Wiechert  (H-W)  formula  for  the  analytical  inversion  of  (2)  and  (3)  is 


-  p/Kr/v) 

In  [ro/r(v)]  =  ^  J  coshT^  (pv/r)  dA 


(4) 


The  inversion  is  applicable  to  cases  in  which  the  derivative  of  A{p)  is  discontinuous,  but  not  to 
cases  in  which  A{p)  itself  is  discontinuous.  Thus  the  Herglotz-Wiechert  formula  cannot  be  used 
in  cases  involving  low-velocity  zone  (LVZ).  However,  Slichter  (1932)  found  an  upper  bound  for 
the  thickness  of  low-velocity  zone  from  observed  gaps  dA  and  ST  in  the  travel  time  curve. 

Once  we  have  the  A  =  A{p)  relationship,  we  can  determine  the  velocity  profile  v  =  v(r) 
and  hence  the  one-dimension  model  can  be  obtained.  It  is  difficult  to  fit  the  slowness  data  shown 
in  Figure  7  using  only  one  curve.  We  separate  the  slowness  data  to  several  branches  which  can  be 
fitted  by  a  set  of  polynomial  curves.  Because  we  do  not  yet  have  not  enough  data  to  determine  the 
local  crust  model,  we  use  the  global  lASPEI  91  cmst  model  as  the  local  model. 


Upper  Mantle  Velocity  Model  for  Western  China 

Unfortunately  because  of  the  scatter  in  the  slowness  measurement,  the  direct  determination 
of  velocity  by  H-W  inversion  of  these  data  alone  does  not  appear  warranted.  Using  the  model 
from  H-W  inversion  as  an  initial  model,  together  with  trial  and  error,  we  determine  a  velocity 
model  which  explains  the  observed  features  and  which  does  not  produce  features  that  are  not 
warranted  by  the  observed  data.  For  example,  the  amplitude  should  provide  important  constraints 
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on  the  velocity  model,  such  as  the  nature  of  the  velocity  discontinuities;  the  ends  of  the  major 
prograde  branches  terminate  with  low  amplitude,  whereas  high  amplitudes  are  associated  with  the 
ends  of  retrograde  branches. 

The  model  W  China  that  we  obtained  is  shown  in  Figure  8  and  is  tabulated  in  Table  2.  The 
model  has  two  well  known  discontinuities  at  417  km  and  689  km.  The  model  also  shows  a  rapid 
increase  in  velocity  between  550  to  600  km  followed  by  low  velocity  gradients  from  600  to  690 
km.  The  travel  times  calculated  from  the  W  China  model  are  compared  with  the  observed 
waveform  data  in  Figure  9. 


Conclusion 

The  use  of  an  array  to  determine  slowness,  and  to  associate  later  arrivals  with  particular 
branches  of  the  P  travel  time  curve  between  10°  and  32°,  has  revealed  a  number  of  consistent 
features.  The  most  striking  feature  of  the  data  is  a  pronounced  later-arrival  branch  from  20°  to  24° 
which  we  associate  with  an  abrupt  velocity  increase  near  a  depth  of  550  km.  Large-amplitude  later 
arrivals  are  also  observed  around  15°,  which  indicate  a  region  of  high  velocity  gradient  near  417 
km.  Both  of  these  regions  of  high  amplitude  later  arrivals  are  on  the  retrograde  ends  of  the 
triplications  resulting  from  velocity  transitions  near  417  km  and  690  km.  Although  determinations 
of  the  exact  depth  and  shapes  of  the  velocity  transitions  are  limited  by  the  non-uniqueness  of  the 
inversion  of  the  scatter  in  our  data,  we  consider  the  main  characteristics  of  the  model  W  China  to 
be  representative  of  the  structure  in  the  upper  mantle  transition  zone  beneath  western  China.  The 
model  shows  a  rapid  increase  in  velocity  between  550  to  600  km  and  low  velocity  gradients  from 
600  to  690  km. 
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Table  1.  Events  used  in  this  study*. 


N  Date  Origin  Time  Lat.  Long,  depth  Dist  Az 

dd/min/yr  hh:rmn:sec  O  (°)  (km)  (°)  O 


1 

02/12/95 

10:56:58.8 

2 

03/05/95 

07:00:00.4 

3 

03/18/95 

18:02:36.6 

4 

04/14/95 

19:53:34.5 

5 

04/24/95 

16:13:11.4 

6 

04/26/95 

03:46:32.8 

7 

04/26/95 

20:11:20.8 

8 

04/26/95 

22:46:27.0 

9 

05/02/95 

11:48:11.6 

10 

05/08/95 

18:37:33.7 

11 

05/12/95 

03:43:49.0 

12 

05/19/95 

11:45:37.3 

13 

05/23/95 

17:17:29.0 

14 

05/30/95 

23:43:51.8 

15 

06/03/95 

11:29:25.2 

16 

06/10/95 

05:08:55.7 

17 

06/11/95 

23:39:51.1 

18 

06/12/95 

03:37:30.0 

19 

06/12/95 

14:47:58.3 

20 

06/21/95 

08:28:41.8 

21 

06/27/95 

05:48:25.1 

22 

06/29/95 

23:03:46.2 

23 

07/04/95 

18:29:45.4 

24 

07/05/95 

23:38:48.9 

25 

07/09/95 

15:56:28.2 

26 

07/09/95 

20:31:31.4 

27 

07/09/95 

20:33:32.5 

28 

07/09/95 

21:29:14.7 

29 

07/10/95 

22:43:16.0 

30 

07/11/95 

21:46:39.7 

31 

07/11/95 

22:20:33.4 

32 

07/11/95 

22:37:07.1 

33 

07/1 1/95 

23:16:16.8 

34 

07/12/95 

03:57:00.4 

35 

07/12/95 

06:02:56.1 

36 

07/12/95 

08:14:33.8 

37 

07/12/95 

09:44:53.1 

38 

07/14/95 

07:29:36.9 

39 

07/16/95 

09:27:12.5 

40 

07/18/95 

11:14:07.2 

41 

07/18/95 

21:33:41.7 

42 

07/19/95 

16:08:22.6 

43 

07/21/95 

22:44:04.5 

44 

07/22/95 

00:06:06.3 

45 

07/23/95 

05:39:49.5 

33.280 

93.386 

27 

24.670 

99.965 

33 

42.422 

87.199 

22 

29.613 

115.616 

10 

22.736 

102.906 

33 

28.895 

103.778 

33 

32.700 

102.388 

33 

37.248 

105.718 

33 

43.776 

84.660 

33 

33.329 

103.592 

33 

23.688 

100.415 

33 

41.720 

92.894 

11 

32.455 

93.560 

33 

27.141 

95.659 

33 

28.080 

102.897 

33 

31.685 

104.208 

33 

31.887 

92.446 

33 

31.781 

104.256 

10 

39.218 

95.292 

24 

37.729 

106.422 

33 

26.209 

96.342 

33 

21.925 

98.943 

33 

30.386 

94.797 

33 

42.513 

86.678 

33 

35.976 

100.073 

33 

21.984 

99.159 

10 

21.725 

99.086 

10 

21.885 

99.116 

10 

21.984 

99.119 

10 

21.966 

99.196 

13 

22.690 

100.099 

10 

21.808 

98.945 

10 

21.467 

98.548 

10 

21.123 

99.034 

10 

21.886 

99.157 

10 

22.041 

99.309 

10 

21.841 

99.161 

10 

21.806 

99.115 

10 

30.317 

94.828 

33 

21.844 

99.038 

10 

22.005 

99.138 

10 

30.321 

94.866 

33 

36.427 

103.123 

13 

36.414 

103.244 

10 

22.015 

98.907 

10 

4.8 

20.48 

332.3 

4.4 

30.76 

333.1 

5.2 

10.12 

327.1 

4.0 

34.72 

317.8 

4.7 

33.72 

331.9 

4.9 

28.85 

325.9 

4.1 

25.06 

322.7 

4.2 

23.50 

313.4 

5.5 

08.01 

330.9 

4.7 

25.19 

320.9 

4.8 

31.82 

333.4 

4.2 

13.31 

316.9 

4.8 

21.28 

333.0 

4.9 

26.86 

335.6 

4.3 

29.10 

327.4 

-4.7 

26.80 

322.4 

4.2 

21.38 

335.3 

4.4 

26.75 

322.3 

5.1 

16.39 

319.6 

4.5 

23.58 

312.1 

4.6 

27.96 

335.6 

4.9 

32.84 

335.9 

4.0 

23.60 

333.6 

4.5 

09.84 

328.3 

5.2 

21.31 

320.1 

5.7 

32.87 

335.7 

5.0 

33.07 

335.9 

4.9 

32.94 

335.8 

4.7 

32.85 

335.7 

6.1 

32.90 

335.6 

4.2 

32.59 

334.3 

4.5 

32.95 

336.0 

4.4 

33.11 

336.6 

4.2 

33.60 

336.3 

4.3 

32.95 

335.7 

4.2 

32.87 

335.5 

4.3 

33.00 

335.8 

4.6 

33.01 

335.8 

4.9 

23.67 

333.7 

4.0 

32.95 

335.9 

4.8 

32.84 

335.7 

4.9 

23.68 

333.6 

5.7 

22.61 

316.6 

4.2 

22.68 

316.5 

4.4 

32.74 

335.9 

11 


Table  1  Continued 


46 

07/25/95 

15:41:24.5 

23.196 

105.845 

10 

4.2 

34.66 

47 

07/30/95 

13:03:59.9 

32.142 

101.440 

33 

4.4 

25.03 

48 

07/30/95 

19:43:09.2 

21.987 

99.052 

33 

4.1 

32.82 

49 

08/02/95 

01:46:31.6 

26.741 

103.053 

33 

4.4 

30.31 

50 

08/02/95 

11:59:43.9 

41.631 

88.447 

10 

4.1 

11.30 

51 

08/07/95 

03:27:16.6 

21.931 

99.218 

33 

4.1 

32.94 

52 

08/17/95 

00:59:57.7 

41.559 

88.800 

0 

6.0 

11.50 

53 

08/22/95 

00:46:53.3 

30.257 

97.203 

33 

4.2 

24.68 

54 

08/28/95 

11:46:05.6 

38.186 

89.001 

33 

4.9 

14.48 

55 

08/30/95 

12:52:59.6 

44.098 

84.117 

33 

4.5 

07.54 

56 

09/02/95 

18:25:26.5 

37.956 

91.395 

33 

4.5 

15.61 

57 

09/04/95 

08:55:35.6 

26.535 

101.638 

18 

4.2 

29.84 

58 

09/04/95 

18:43:45.1 

43.903 

87.436 

33 

4.1 

09.02 

59 

09/05/95 

16:19:17.0 

27.094 

96.680 

33 

4.3 

27.28 

60 

09/06/95 

00:11:23.1 

29.156 

104.012 

33 

4.3  . 

28.75 

61 

10/06/95 

14:26:04.7 

33.407 

93.467 

33 

4.4 

20.40 

62 

10/09/95 

17:48:23.8 

25.075 

97.727 

33 

4.4 

29.51 

63 

10/14/95 

04:50:32.5 

33.158 

92.783 

33 

4.2 

20.36 

64 

10/14/95 

05:52:46.5 

44.422 

86.926 

33 

4.0 

08.39 

65 

10/21/95 

20:10:49.0 

36.215 

103.301 

33 

4.4 

22.86 

66 

10/23/95 

22:46:50.8 

26.003 

102.227 

10 

5.8 

30.57 

67 

10/23/95 

23:47:10.3 

25.791 

102.561 

33 

4.5 

30.90 

68 

10/24/95 

00:09:54.7 

25.846 

102.146 

10 

4.2 

30.67 

69 

10/24/95 

00:29:03.5 

25.847 

102.106 

10 

4.3 

30.65 

70 

10/24/95 

01:00:44.0 

25.858 

102.223 

10 

4.3 

30.69 

71 

10/24/95 

01:19:11.5 

25.830 

102.236 

10 

4.3 

30.72 

72 

10/24/95 

01:40:01.7 

25.923 

102.365 

10 

4.2 

30.70 

73 

10/24/95 

02:11:59.4 

25.801 

102.067 

10 

4.3 

30.67 

74 

10/24/95 

08:08:14.4 

25.886 

102.547 

10 

4.6 

30.81 

75 

10/24/95 

10:45:07.9 

25.719 

102.281 

33 

4.2 

30.84 

76 

10/24/95 

17:29:55.9 

25.754 

102.098 

10 

4.4 

30.73 

77 

10/24/95 

17:43:26.1 

25.914 

102.340 

10 

4.7 

30.70 

78 

10/24/95 

22:32:54.6 

25.906 

102.163 

10 

4.7 

30.62 

79 

10/24/95 

23:13:15.8 

25.777 

102.110 

10 

4.4 

30.71 

80 

10/25/95 

02:42:57.0 

25.903 

102.430 

10 

4.3 

30.75 

81 

10/25/95 

03:16:32.3 

25.963 

102.210 

10 

4.4 

30.59 

82 

10/25/95 

04:21:52.1 

25.773 

102.256 

10 

4.4 

30.78 

83 

10/25/95 

17:33:39.0 

25.879 

102.189 

10 

4.1 

30.66 

84 

10/26/95 

13:19:34.9 

25.738 

102.045 

10 

4.0 

30.72 

85 

10/27/95 

02:05:36.4 

25.896 

102.201 

10 

4.3 

30.65 

86 

10/27/95 

08:52:34.6 

36.696 

92.080 

33 

4.9 

16.97 

87 

10/28/95 

12:55:05.2 

26.020 

102.300 

10 

4.4 

30.59 

88 

10/28/95 

20:36:03.7 

25.958 

102.257 

10 

4.5 

30.62 

89 

10/29/95 

15:46:01.5 

25.667 

102.083 

33 

4.2 

30.80 

90 

10/29/95 

17:18:42.3 

25.852 

102.219 

10 

4.5 

30.70 

91 

11/01/95 

08:10:07.8 

25.868 

102.339 

10 

4.8 

30.74 

92 

11/01/95 

15:29:36.2 

25.926 

102.247 

11 

4.5 

30.64 

93 

11/02/95 

21:26:51.2 

25.836 

102.069 

10 

4.0 

30.64 

94 

11/04/95 

20:07:39.8 

33.018 

92.140 

33 

4.0 

20.24 

95 

11/05/95 

17:14:53.6 

32.896 

92.201 

24 

5.1 

20.38 

96 

11/07/95 

03:50:09.4 

26.004 

102.471 

10 

4.3 

30.68 

97 

11/07/95 

20:47:07.2 

25.872 

102.425 

20 

4.6 

30.77 

329.2 

324.3 

335.8 

328.6 

326.1 

335.7 

325.4 

330.9 

332.5 

331.9 

328.3 
330.0 

321.2 

334.5 

325.4 
332.0 

334.9 

333.3 

320.6 

316.8 

329.9 

329.8 

330.1 

330.1 
330.0 

330.1 

329.9 

330.2 

329.7 

330.1 

330.2 

329.9 

330.1 

330.2 

329.8 
330.0 

330.1 

330.1 

330.3 
330.0 

329.4 

329.8 

329.9 

330.3 

330.1 

329.9 
330.0 

330.2 

334.4 

334.5 
329.7 

329.9 


12 


Table  1  Continued 

98  11/09/95  20:55:17.6 

22.058 

102.724 

33 

4.4 

34.24 

332.5 

99 

11/11/95 

14:29:27.0 

34.935 

104.080 

10 

4.3 

24.23 

318.2 

100 

11/18/95 

04:48:11.7 

25.673 

102.342 

20 

4.0 

30.91 

330.1 

101 

11/29/95 

14:23:08.7 

25.425 

102.142 

10 

4.0 

31.03 

330.5 

102 

12/03/95 

16:17:23.8 

25.984 

102.193 

20 

4.1 

30.57 

330.0 

103 

12/07/95 

01:28:50.1 

36.406 

100.114 

33 

4.0 

21.00 

319.4 

104 

12/12/95 

05:42:32.2 

26.852 

96.059 

33 

5.1 

27.27 

335.4 

105 

12/12/95 

17:31:16.8 

42.117 

86.911 

33 

4.3 

10.27 

328.8 

106 

12/18/95 

04:57:04.7 

34.511 

97.351 

33 

5.3 

21.13 

325.4 

107 

12/20/95 

10:07:58.7 

34.488 

97.474 

33 

5.0 

21.20 

325.3 

108 

12/26/95 

03:31:09.7 

29.313 

105.453 

33 

4.2 

29.35 

324.1 

109 

01/02/96 

11:37:41.1 

38.515 

102.845 

33 

4.2 

20.98 

313.1 

no 

01/04/96 

02:55:45.5 

38.140 

102.815 

33 

4.1 

21.22 

313.9 

111 

01/04/96 

03:49:27.6 

38.740 

104.598 

33 

4.7 

21.84 

311.5 

112 

01/04/96 

07:26:52.4 

36.218 

100.491 

33 

4.2 

21.34 

319.3 

113 

01/09/96 

06:27:54.4 

43.703 

85.654 

33 

5.2 

08.44 

327.6 

114 

01/16/96 

00:08:18.5 

29.458 

105.605 

33 

4.2 

29.31 

323.9 

115 

01/18/96 

00:19:54.5 

23.894 

94.597 

33 

4.0 

29.48 

339.1 

116 

01/26/96 

02:21:11.2 

30.873 

91.478 

33 

5.1 

21.98 

337.7 

117 

02/02/96 

22:59:55.5 

22.141 

103.164 

10 

4.2 

34.36 

332.1 

118 

02/03/96 

11:14:19.8 

27.299 

100.341 

10 

6.3 

28.61 

330.5 

119 

02/03/96 

11:24:01.9 

27.280 

100.820 

10 

4.5 

28.83 

330.1 

120 

02/03/96 

11:39:41.9 

27.267 

100.318 

10 

4.9 

28.62 

330.6 

121 

02/03/96 

11:58:24.5 

27.151 

100.331 

10 

4.6 

28.73 

330.6 

122 

02/03/96 

12:03:28.8 

27.431 

100.877 

10 

4.2 

28.73 

329.9 

123 

02/03/96 

12:08:16.6 

27.207 

100.358 

10 

5.0 

28.69 

330.6 

124 

02/03/96 

12:14:52.1 

27.053 

100.465 

10 

4.6 

28.88 

330.6 

125 

02/03/96 

12:36:59.6 

27.165 

100.539 

10 

4.0 

28.81 

330.4 

126 

02/03/96 

12:49:11.1 

27.115 

100.584 

10 

4.0 

28.87 

330.4 

127 

02/03/96 

13:17:30.0 

27.214 

100.254 

10 

4.3 

28.64 

330.7 

128 

02/03/96 

13:35:06.4 

27.150 

100.365 

10 

4.0 

28.75 

330.6 

129 

02/03/96 

14:00:14.7 

27.316 

100.654 

10 

4.2 

28.73 

330.2 

130 

02/03/96 

14:19:36.6 

27.131 

100.379 

10 

4.4 

28.77 

330.6 

131 

02/03/96 

14:56:27.2 

27.372 

100.380 

10 

4.1 

28.56 

330.4 

132 

02/03/96 

15:10:02.5 

26.798 

100.124 

10 

4.0 

28.95 

331.2 

133 

02/03/96 

17:49:39.1 

27.310 

100.724 

10 

4.0 

28.77 

330.1 

134 

02/03/96 

19:38:58.1 

27.079 

100.527 

10 

4.7 

28.88 

330.5 

135 

02/04/96 

04:37:47.9 

26.976 

100.447 

10 

4.1 

28.93 

330.7 

136 

02/04/96 

16:58:05.8 

27.050 

100.383 

10 

5.5 

28.84 

330.7 

137 

02/04/96 

22:57:48.4 

26.972 

100.343 

10 

4.6 

28.89 

330.8 

138 

02/05/96 

03:26:10.5 

27.010 

100.424 

10 

4.7 

28.89 

330.7 

139 

02/05/96 

14:36:06.4 

21.366 

97.724 

33 

4.1 

32.90 

337.4 

140 

02/05/96 

19:40:21.6 

26.898 

100.305 

10 

4.3 

28.94 

330.9 

141 

02/06/96 

01:42:32.5 

27.049 

100.314 

10 

4.3 

28.81 

330.8 

142 

02/06/96 

07:36:13.8 

27.247 

100.418 

10 

5.3 

28.68 

330.5 

143 

02/06/96 

20:08:57.4 

27.241 

100.530 

10 

5.2 

28.74 

330.4 

144 

02/07/96 

07:19:41.4 

26.969 

100.505 

10 

4.8 

28.96 

330.6 

145 

02/07/96 

15:29:23.9 

27.125 

100.894 

10 

4.2 

29.00 

330.1 

146 

02/07/96 

21:04:27.4 

27.247 

100.691 

10 

4.3 

28.80 

330.2 

147 

02/08/96 

15:54:07.9 

25.408 

96.991 

33 

4.2 

28.93 

335.5 

148 

02/08/96 

22:46:40.2 

27.019 

100.499 

10 

4.4 

28.92 

330.6 

13 


Table  1  Continued 

149  02/11/96  08:43:28.2 

27.529 

101.301 

10 

4.4 

28.83 

329.4 

150 

02/15/96 

01:09:50.7 

27.252 

100.606 

33 

4.3 

28.76 

330.3 

151 

02/16/96 

03:02:22.5 

28.863 

92.900 

11 

4.4 

24.30 

337.5 

152 

02/24/96 

00:57:39.0 

27.208 

100.278 

33 

4.3 

28.66 

330.6 

153 

02/28/96 

11:22:02.0 

29.099 

104.737 

33 

5.1 

29.16 

324.9 

154 

02/28/96 

23:21:48.4 

27.045 

100.279 

33 

4.1 

28.80 

330.8 

155 

03/13/96 

03:52:14.3 

27.039 

100.389 

33 

4.7 

28.85 

330.7 

156 

03/16/96 

06:46:53.5 

27.139 

100.625 

33 

4.1 

28.87 

330.4 

157 

03/19/96 

10:48:38.8 

24.946 

96.418 

33 

4.4 

29.14 

336.4 

158 

03/20/96 

02:11:21.9 

42.182 

87.627 

25 

4.8 

10.50 

326.6 

159 

03/21/96 

22:34:04.8 

26.351 

98.450 

10 

4.3 

28.65 

333.2 

160 

03/31/96 

03:07:14.0 

43.019 

88.680 

33 

4.2 

10.28 

321.2 

161 

04/09/96 

22:00:46.1 

20.282 

98.719 

33 

4.0 

34.26 

337.1 

162 

04/10/96 

04:13:18.4 

24.026 

96.714 

33 

4.1 

30.09 

336.7 

163 

04/16/96 

21:24:22.6 

26.659 

100.152 

10 

4.2 

29.08 

331.2 

164 

04/22/96 

04:43:21.6 

30.048 

92.074 

33 

4.5 

22.93 

337.6 

165 

04/22/96 

20:24:40.1 

44.141 

86.925 

33 

4.0 

08.61 

321.8 

Table  2.  W  China  P  velocity  model*. 


Depth 

(km) 

Vp 

(km/s) 

Depth 

(km) 

Vp 

(km/s) 

Depth 

(km) 

Vp 

(km/s) 

0 

5.80 

260 

8.54 

530 

9.65 

10 

5.80 

270 

8.57 

540 

9.70 

20 

5.80 

280 

8.61 

550 

9.76 

20 

6.50 

290 

8.64 

560 

9.84 

30 

6.50 

300 

8.68 

570 

9.96 

35 

6.50 

310 

8.71 

580 

10.10 

35 

8.00 

320 

8.75 

590 

10.21 

40 

8.04 

330 

8.78 

600 

10.32 

50 

8.04 

340 

8.82 

610 

10.34 

60 

8.04 

350 

8.85 

620 

10.36 

70 

8.04 

360 

8.89 

630 

10.37 

80 

8.04 

370 

8.92 

640 

10.38 

90 

8.04 

380 

8.96 

650 

10.40 

100 

8.04 

390 

8.99 

660 

10.40 

110 

8.05 

400 

9.03 

670 

10.42 

120 

8.06 

410 

9.07 

680 

10.45 

130 

8.09 

417 

9.13 

689 

10.55 

140 

8.12 

417 

9.30 

689 

10.90 

150 

8.16 

420 

9.31 

690 

10.92 

160 

8.19 

430 

9.34 

700 

10.96 

170 

8.22 

440 

9.37 

710 

10.99 

180 

8.26 

450 

9.41 

720 

11.03 

190 

8.29 

460 

9.45 

730 

11.06 

200 

8.33 

470 

9.49 

740 

11.10 

210 

8.36 

480 

9.52 

750 

11.14 

220 

8.40 

490 

9.56 

760 

11.16 

230 

8.43 

500 

9.58 

770 

11.17 

240 

8.47 

510 

9.60 

780 

11.18 

250 

8.50 

520 

9.62 

790 

11.19 

*  The  crust  and  upper  most  mantle  model  is  lASPEI  91  model. 
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Figure  Captions 


Fig.  1.  Location  of  the  cross-array  in  Kurchatov,  Kazakstan.  The  configuration  of  the  21 -element 
array  with  a  station  spacing  of  2.25  km  is  depicted.  Other  seismographic  stations  in  central  Asia 
are  also  shown  for  reference. 

Fig.  2.  Epicenters  of  the  earthquakes  used  in  this  study  are  plotted  on  a  map  of  topographic  relief. 
The  size  of  the  circle  represents  the  size  of  the  earthquake,  and  the  focal  mechanisms  (Harvard 
CMT  solutions)  of  the  earthquakes  in  the  region  are  plotted  as  beach  balls.  P-wave  signals  from 
most  of  the  earthquakes  selected  for  analysis  propagate  through  the  crust  and  upper  mantle  beneath 
the  Tibetan  Plateau  andthe  Northern  Tien-Shan  mountains  in  western  China. 

Fig.  3.  Seismic  record  section  showing  the  waveforms  from  earthquakes  in  the  distance  range  10° 
to  20°.  Records  at  the  center  of  the  array  (C21)  are  plotted  with  a  reduction  velocity  of  10  km/s 
(=11.12  s/deg)  and  predicted  travel  time  curves  for  focal  depth  =  15  km  for  various  P-wave 
branches  are  indicated.  The  A-B  branch  is  for  Pn  arrivals,  the  C-D  branch  for  first  arriving  P 
waves  whose  rays  have  turning  point  below  the  410  km  discontinuity,  and  G-H  is  a  portion  of  the 
P  travel  time  branch  corresponding  to  rays  turning  below  the  "660  km  discontinuity"  (later  in  this 
paper,  estimated  as  having  a  depth  of  about  680  km).  Note  that  record  from  each  event  is  plotted 
here  without  adjusting  for  its  focal  depth. 

Fig.  4.  A  synthetic  record  section  showing  the  first  as  well  as  later  arrival  P  waves  in  the  lASPEI 
91  model  for  a  focal  depth  of  15  km.  The  source  mechanism  from  an  event  on  12/18/95  (04:57, 
M=5.7)  is  used  to  calculate  amplitude  and  travel  times  using  the  WKBJ  seismogram  method. 

Fig.  5.  An  example  of  the  method  of  determination  of  slowness  and  azimuth  by  using  a  weighted, 
delay-and-sum  beamform.  The  event  analyzed  occurred  on  05/23/95  (17:17,  A  =  21.28°,  /n^  = 
4.8,  h  =33  km,  PDE).  (a)  the  top  trace  is  the  reference  station  waveform  and  the  bottom  trace  is 
the  array  beam  calculated  using  the  slowness  and  azimuth  determined  from  the  time  window 
indicated;  (b)  azimuth  (cross),  slowness  (circle)  and  beam  power  (triangle)  determined  at  each 
time  window  are  plotted  vs  time;  (c)  contour  of  the  beam  power  for  the  best  beam  determination; 
(d)  waveforms  at  all  stations  are  plotted,  each  with  its  appropriate  delay,  and  the  top  trace  is  the 
array  beam. 

Fig.  6.  (a)  Combined  first  arrival  travel  time  and  slowness  data  reduced  to  a  surface  focus  are 
plotted  with  a  reduction  slowness  of  11.12  s/deg  (thick  lines),  and  showing  corresponding 
slownesses  also  corrected  for  the  local  structure  (thin  lines).  Slope  lines  plotted  using  the  slope  of 
the  measured  slowness,  (b)  The  slownesses  of  the  first  arrivals  corrected  to  surface  focus  are 
plotted  vs.  distance  (open  circles),  and  showing  corresponding  slownesses  corrected  also  for  the 
local  structure  (solid  circles). 

Fig.  7.  Slownesses  corrected  for  local  structure  and  reduced  to  the  surface  focus  are  plotted  against 
epicentral  distance.  The  slownesses  are  fitted  by  several  polynomial  curves  and  travel  time 
branches  are  marked  by  letters. 
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Fig.  8.  (a)  Slownesses  corrected  for  local  structure  and  reduced  to  the  surface  focus  are  plotted 
against  epicentral  distance.  The  slownesses  are  fitted  by  several  polynomial  curves  identified  by 
corresponding  P-wave  travel  time  branches,  (b)  W  China  upper  mantle  model  (solid  line) 
obtained  in  this  study  is  compared  with  the  lASPEI  91  model  (dashed  line). 

Fig.  9.  Selected  waveform  data  are  plotted  with  reduced  velocity  of  11.12  s/deg  and  travel  time 
curves  predicted  by  the  W  China  model  are  superimposed  for  comparison. 
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